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摘要 : 针对 小 波 变 换 电 压 暂 降 检 测算 法 在 谐 波 干扰 条 件 下 的 不 足 ， 提 出 了 一 种 基 
于 背 窗 小 波 和 二 点 法 的 电压 暂 降 算法 。 该 算法 首先 通过 滑 窗 小 波 滤 除 暂 降 信 号 中 的 谐 
波 分 量 ， 得 到 电网 基 波 信号 ， 再 利用 二 点 法 计算 电压 暂 降幅 值 比 及 相位 跳 变 角 。 实 例 
仿真 表明 ， 所 提 算 法 不 存在 因 小 波 系数 最 大 模 难 以 识别 而 导致 的 暂 降 特 征 量 无 法 计算 
的 问题 ， 能 够 在 谐 波 干扰 条 件 下 准确 地 检测 出 电压 暂 降 特 征 参数 ， 且 相关 检测 结果 延 
述 小 于 一 个 工 频 周 期 ,满足 实际 工程 需 
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Study on Detection of Voltage Sag Based on Sliding Window 
Wavelet and Two Points Method 
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( Nanjing University of Science and Technology Nanjing 210094 China) 


Abstract: A algorithm of voltage sag detection based on sliding window wavelet 
and two points method was proposed for overcome the disadvantage of wavelet transform 
in the condition of harmonics interference. According to the proposed algorithm, the 
fundamental voltage signal can be obtained by using sliding window wavelet to filter 
harmonic component. Then, the amplitude ratio and phase jump information of voltage 
sag are detected by Two Points method. Simulation results show that the problem of 
voltage sag characteristic parameter calculation, which caused by wavelet coefficient 
maximum norm is hardly to identification, does not exist when the proposed algorithm is 
adopted. The proposed algorithm can meet the needs in engineering due to the voltage sag 
characteristic parameter can be detected under the condition of harmonics interference and 
the delay of detecting results less than one cycle. 
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1 引言 


随 着 电力 负荷 增长 ， 电 压 暂 降 引 起 的 电网 暂 态 
扰动 已 引起 广泛 关注 。 尤 其 在 配 网 中 ， 电 压 暂 降 
导致 的 电网 故障 可 占 故 障 总 数 的 半数 以 上 。 因 此 ， 
有 必要 研究 准确 的 电压 暂 降 检测 算法 ， 以 实现 暂 降 
源 的 正确 判断 和 电能 质量 的 有 效 改善 。 

目前 ， 电 压 暂 降 检测 的 常用 算法 有 RMS 法 、 
求 导 法 、 瞬 时 无 功率 法 和 短 时 傅 里 叶 变 换 法 等 。 其 
中 ，RMS 法 通过 计算 每 个 周波 电压 有 效 值 来 判断 暂 
降 的 发 生 ， 是 一 种 经 典 检测 法 ， 有 具有 原理 简单 ， 计 
算 量 小 等 优势 。 但 该 方法 对 谐 波 和 噪声 较为 敏感 ， 
精度 有 限 ， 且 无 法 描述 暂 降 后 的 电压 相 角 变化 "。 
求 导 算 法 同样 具有 计算 量 小 、 易 实现 的 优势 ， 但 其 
计算 精度 亦 会 受到 谐 波 和 噪声 干扰 ， 且 还 会 受到 采 
样 频率 的 影响 “”"。 瞬 时 无 功 功率 法 在 dq 坐标 或 aB 
坐标 中 计算 电压 暂 降 参数 ， 精 度 较 高 ， 但 在 单 相 或 
三 相 不 对 称 系统 中 应 用 时 存在 一 定 相 位 延迟 ， 且 算 
法 需 进行 低 通 滤 波 ， 实 时 性 会 进一步 下 降 "。 短 时 
傅 里 叶 变 换 法 将 非 平稳 信号 视 为 一 系列 短 时 平稳 信 
号 的 全 加 ， 其 短 时 性 通过 时 间 窗 来 获得 ， 在 合适 的 窗 
函数 条 件 下 ， 所 得 结果 具有 较 高 精度 和 实时 性 """。 该 
算法 缺陷 在 于 窗口 形状 一 经 确定 便 无 法 改变 ， 这 相 
当 于 使 用 固定 尺度 对 信号 进行 分 析 ， 故 对 于 具有 有 瞬 
变 特性 的 电压 暂 降 信号 ， 窗 函数 选择 较为 困难 一 ， 
具有 一 定局 限 性 。 

近年 来 ， 小 波 变 换 方法 被 广泛 应 用 于 电能 质量 
分 析 问题 中 一 。 小 波 变换 可 通过 尺度 参数 和 位 移 
参数 的 共同 变化 ， 改 变 小 波 窗口 的 大 小 及 形状 。 信 
号 频带 越 高 ， 分 析 窗 口 越 小 ， 时 域 分 辩 率 越 高 。 因 
此 ， 小 波 变 换 法 的 优势 在 于 能 在 时 频 两 域 良 好 反映 
信号 的 局 部 特性 ， 准 确 提取 信号 特征 ， 尤 其 适合 分 
析 电 网 的 暂 态 扰动 。 

然而 ， 小 波 变 换 在 电压 暂 降 检测 中 应 用 时 ， 也 
存在 一 定 缺陷 。 这 是 因为 电压 暂 降 时 刻 通 常 利用 小 
波 分 解 后 的 高 频 小 波 系数 或 其 模 极 大 值 来 确定 ， 电 
压 暂 降幅 值 和 相位 等 特征 参量 也 由 此 得 到 '"。 但 
是 ， 电 网 电压 信号 往往 含有 谐 波 分 量 ， 而 电压 暂 降 
发 生 时 刻 亦 是 随机 的 ， 故 带 有 谐 波 的 电压 暂 降 信号 
频率 特性 较为 复杂 。 此 时 ， 电 压 暂 隆 引 起 的 小 波 系 
数 变 化 可 能 难以 辨别 ， 造 成 暂 降 定 位 能 力 减弱 "， 
进而 导致 暂 降 特征 量 计算 困难 。 因 此 ， 有 必要 对 基 
于 小 波 变 换 的 电压 暂 降 检测 法 加 以 改进 ， 来 克服 上 
述 缺 陷 。 
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为 此 ， 本 文 提 出 了 基于 请 窗 小 波 和 二 点 算法 的 
电压 暂 降 检测 方法 ， 该 方法 并 不 利用 高 频 小 波 系数 
检测 暂 降 参 数 ， 而 是 利用 请 窗 小 波 对 电压 信号 进行 
滤波 ， 排 除 谐 波 干扰 ， 再 通过 二 点 算法 计算 电压 暂 
降 比 及 相 角 变化 量 。 仿 真 结果 表明 ， 该 方法 能 够 在 
谐 波 干扰 条 件 下 有 效 检测 电压 暂 降 过 程 。 


2 ”基于 小 波 变换 的 电压 暂 降 检测 

小 波 变换 可 将 信号 () 通过 不 同 尺 度 分 解 至 尺 
度 空间 和 小 波 空间 中 。 此 时 ， 信 号 可 用 多 分 辩 率 公 
式 表示 

JDO=-Zow(Dex(D+ 之 2d0DWxO (0 


式 中 ，a( 有 为 概貌 系数 ，d/( 忆 为 小 波 系 数 ， PD 
为 尺度 函数 ， 到 AD 为 小 波 函 数 ， 下 标 j 为 小 波 分 解 
层 数 。 

式 (1) 中 ，wxD 和 到 AD) 为 尺度 空间 和 小 波 
空间 的 标准 正 交 基 。a(D 和 4d(D 可 由 如 下 分 解 公式 
求解 


aja(k)= Hz-26ai(D 
dmua(D = 之 SO-2D4i( 
式 中 ，h(n) 和 g(n) 分 别 为 小 波 变换 的 低 通 和 高 通 滤 
波 器 系数 。 
小 波 分 解 后 ， 第 7 层 的 概貌 系数 可 由 第 +1 层 
的 概貌 系数 和 小 波 系数 重 构 得 到 


aj(m)= hn-2pan(k)+ Seln-20)dja(k) (2) 


上 述 小 波 变换 分 解 及 重 构 过 程 的 示意 图 如 图 1 
所 示 。 


0 


1 ”小 波 分 解 及 重 构 示意 图 
Fig.l Schematic diagram of wavelet decomposition and 


wavelet reconstruction 


经 过 合理 的 层 数 划分 和 采样 频率 选取 ， 小 波 变 
换 的 各 层 小 波 系数 能 够 有 效 反 映 电网 电压 谐 波 分 量 。 
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当 电 网 电压 发 生 暂 降 时 ， 原 始 电压 信号 中 等 效 于 
全 加 了 一 个 高 频 分 量 ， 因 此 通过 检测 小 波 系 数 的 突 
变 情况 可 以 有 效 定位 暂 降 时刻， 从 而 实现 电压 暂 降 
检测 。 

然而 ， 电 网 电压 暂 降 时 刻 是 随机 的 ， 电 网 谐 波 
分 量 的 频率 与 相位 亦 是 随机 出 现 的 ， 因 此 暂 降 引起 
的 高 频 分 量 可 能 与 谐 波 分 量 相 互 干 扰 ， 造 成 小 波 系 
数 突变 点 难以 识别 ， 进 而 导致 电压 暂 降 检测 精度 下 
降 ， 甚 至 检测 失效 。 此 为 小 波 变换 电压 暂 降 检测 的 
主要 缺陷 。 
3 ”基于 滑 窗 小 波 和 二 点 算法 的 电压 暂 降 

检测 
3.1 电压 暂 降 检测 二 点 法 原理 

为 克服 小 波 变 换 电压 暂 降 检测 的 不 足 ， 首 先 提 


出 一 种 暂 降 检测 的 二 点 法 ， 下 面 推导 算法 原理 。 
设 单 相 电压 暂 降 前 为 理想 电压 u(t)， 则 有 


u(t)= V2U sin (wt+9) 
式 中 ,UU 为 电压 有 效 值 ，w 为 角 频 率 ， 9 为 初 相 角 。 
暂 降 后 电压 变 为 
u(t)= kV2U sin (wt+9+09) 
式 中 ,为 电压 暂 降 的 幅 值 比 ; 少 为 电压 暂 降 的 相位 
跳 变 角 。Kk 和 9 即 为 电压 暂 降 特征 量 。 


设 相 邻 两 个 采样 时 刻 的 电网 电压 信号 za(D) 及 
za(D 表达 式 为 


u(t)= kV2U sin (wt+m+9) 
w(t) = KV2U sin (wt +A)+O+9) 


其 中 
今 
x -ee (3) 
Xx, ksing 
再 令 
CQ =CcoS (OA1) 
P=sin (WAt) 


A= V2sin (wt +9) 
B= V2cos (wt+9) 


式 中 ,4 和 B 可 由 锁 相 环 得 到 。 
将 wi(?) 和 w(?) 进行 三 角 变 换 ， 得 
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HX=U (4) 


oo 
w(t)/U 


4 B 
CD 
C= Aa+BPp 
D=Ba- 4B 
压 暂 降 特 征 量 的 求解 。 
下 面 分 析 一 下 式 (4) 的 有 解 条 件 。 若 系数 和 矩阵 
五 不满 秩 ， 即 


U= 


AD-BC=0 


由 式 (4) 可 得 增 广 和 矩阵 


4 B CD/U 


[HIUV]= 
4a Ba w(t)/U 


(5) 


根据 a 和 的 定义 ， 当 B=0 时 ， 必 有 Q=1。 
那么 ， 当 m(D=z(0 时 ， 增 广 矩 阵 与 系数 阵 等 秩 ， 
但 不 满 秋 ， 式 (4) 有 无 穷 解 ， 当 m(O <=z(0 时 ， 增 
广 和 矩阵 与 系数 阵 不 等 秋 ， 式 (4) 无 解 。 

当 且 仅 当 Bz=0 时 ， 式 (5) 有 唯一 解 


x __Bu (Da- Au()B- Bu,(t) 
20B 
1 AD- Au Da- Bu()B 
- 20B 


根据 式 (3)， 可 得 电压 暂 降 特征 量 


k=VX + 
| (6) 


Xl 


0= | 


根据 Bz0 条 件 ， 易 得 
At zn/(27) (7) 


式 中 , f 为 电网 频率 ,n=0,1,2…。 可 见 ， 只 要 采样 
周期 满足 式 (7)， 即 可 求 得 特征 量 。 对 于 电网 工 频 
信号 而 言 ， At 通常 远 小 于 f， 不 构成 无 穷 解 条 件 ， 
而 满足 唯一 解 条 件 。 

从 算法 原理 上 看 ， 二 点 法 仅 需 两 个 采样 点 即 可 
计算 电压 暂 降 特征 量 ， 具 有 良好 快速 性 。 但 由 于 算 
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法 是 基于 理想 正弦 电压 推导 的 ， 因 此 电网 谐 波 必然 
会 降低 其 计算 精度 。 


3.2 滑 窗 小 波 变 换 与 二 点 算法 结合 的 电压 暂 降 检测 
原理 


谐 波 会 对 电压 暂 降 检测 的 小 波 变换 法 和 二 点 法 
的 精度 造成 影响 ， 为 此 提出 一 种 基于 清 窗 小 波 和 二 
点 法 的 电压 暂 降 检测 算法 。 

该 算法 的 思想 是 利用 小 波 变 换 的 时 频 域 分 析 能 
力 滤 除 电 压 谐 波 干扰 ， 提 取 基 波 信号 ， 再 通过 二 点 
法 对 暂 降 特征 量 进行 检测 。 由 于 小 波 变 换 需要 一 个 
工 频 周期 采样 数据 参与 运算 ， 这 意味 着 电压 暂 降 信 
号 的 基 波 检测 将 至 少 延迟 一 个 工 频 周期 ， 因 此 可 在 
小 波 变换 中 加 入 滑动 时 间 窗 ， 使 之 构成 请 窗 小 波 “"。 
请 窗 小 波 能 在 时 域 中 对 采样 信号 进行 滑动 加 窗 ， 使 
每 个 时 刻 的 电压 采样 值 均 与 之 前 采样 值 构成 小 波 变 
换 初 始 信号 ， 从 而 实现 电压 暂 降 特征 量 的 检测 延迟 
小 于 一 个 工 频 周期 。 

基于 清 窗 小 波 和 二 点 法 的 电压 暂 降 检测 实现 流 
程 如 图 2 所 示 。 


开始 


Y 


获得 第 7 时 刻 电 压 采样 值 | 


Y 
利用 滑动 窗 将 第 i 时 刻 采样 值 与 之 前 时 刻 
的 采样 值 构造 小 波 变 换 初始 信号 


Y 
进行 小 波 变换 ， 得 到 基 波 信号 | 


Y 
任意 选取 基 波 信号 中 的 两 个 时 刻 值 进行 
二 点 法 电压 暂 降 检测 


得 到 第 ;时 刻 的 电压 暂 降 的 幅 值 比 大 及 
tk 变 4 


Wal 


图 2 基于 滑 窗 小 波 和 二 点 法 的 电压 暂 降 检测 实现 流程 
Fig.2 Implementation procedure for voltage sag detection 


based on sliding window wavelet and two points method 


4 ”实例 仿真 与 结果 分 析 


4.1 仿真 信号 
设 一 组 带 有 谐 波 的 电压 暂 降 仿真 信号 如 下 
u(t) = M2U,, sin (Ot+ pi) +30V2U,, sin (ot -nx/3)+ 
15V2U, sin (Swt +n/6) (8) 
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式 中 ， 为 基 波 电压 幅 值 ，9, 为 基 波 电压 相 角 ，w 
为 角 频率 。 车 电压 暂 降 发 生 时 刻 为 4， 电 压 恢复 时 
刻 为 4，Us 初始 值 为 20V， 9 初始 值 为 -23 rad， 
电压 暂 降 幅 值 比 x 为 0.5， 暂 降 相 位 跳 变 角 9 为 
-7/6 rad， 则 暂 降 发 生前 后 的 电压 幅 值 及 相位 可 表示 
如 下 


220 ft<t,t>t -TV/3 ft<t,t>t 
本 ? Pn 
110 <t<t -xX/2 ti<t<t 


m 


4.2 小 波 变换 电压 暂 降 检测 的 仿真 分 析 

首先 利用 式 (8) 所 示 的 仿真 信号 验证 小 波 变 
换 电压 暂 降 检测 的 不 足 。 为 使 谐 波 和 电压 暂 降 对 念 
真 信号 的 时 频 特 性 造成 影响 ， 仿 真 在 如 表 1 所 示 
的 不 同和 条 件 下 进行 。 仿 真 时 ， 采 样 频率 为 
3.2kHz。 


表 1 电压 暂 降 和 恢复 时 刻 的 取 值 


Tab.1 The time of voltage sag and recovery 


序号 暂 降 时 刻 is 恢复 时 刻 ws 
Casel 0.0766 0.2263 
Case2 0.0781 0.2278 
Case3 0.0797 0.2294 
Case 4 0.0813 0.2309 
Case 5 0.0825 0.2322 


第 一 层 小 波 变 换 后 的 系数 Di 高 频 分 量 最 为 丰 
富 ， 因 此 通常 用 该 层 小 波 系数 最 大 模 值 来 定位 电压 
暂 降 过 程 。 图 3 给 出 了 不 同和 条 件 下 ， 利 用 
dB30 小 波 得 到 的 小 波 系 数 Di 仿真 结果 。 易 见 ， 在 
3a ~ 3d 中， 小 波 系数 在 电压 暂 降 和 恢复 时 刻 均 
有 明显 突变 。 但 在 图 3e 中 ， 小 波 系数 突变 点 却 在 暂 
降 信 号 频率 分 量 与 谐 波 频率 分 量 的 相互 影响 下 ， 变 
得 难以 辨识 ， 因 此 无 法 准确 定位 电压 暂 降 和 恢复 时 
刻 ， 进 而 导致 暂 降 特征 参数 无 法 求解 。 这 验证 了 传 


$50 r 
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(a) t=0.076 6s, 1.=0.226 3s 
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(d) £4=0.081 3s, £.=0.230 9s 


0 00 010 015 020 025 050 
(e) £1=0.082 5s, £.=0.2322s 


3 不 同和 .条件 下 的 小 波 变换 电压 暂 降 检 测 结 果 
Fig.3 The results of voltage sag detection based on wavelet 


transform with the various ty and 大 


统 小 波 变 换 检 测 方 法 在 谐 波 干扰 条 件 下 ， 可 能 存在 
失效 现象 。 
4.3 二 点 法 电压 暂 降 检测 的 仿真 分 析 

再 对 电压 暂 降 检测 二 点 法 的 特点 进行 仿真 分 析 。 
首先 去 除 暂 降 仿真 信号 中 的 谐 波 干扰 ， 在 无 谐 波 条 
件 下 分 析 算 法 优势 。 图 4 中 的 仿真 结果 表明 当 电 压 
发 生 暂 降 和 恢复 时 ， 二 点 法 所 得 暂 降幅 值 比 及 相位 
跳 变 角 检 测 相 对 误差 分 别 仅 为 0.28% 和 0.35%， 且 
几乎 无 延迟 。 这 验证 了 算法 在 电网 无 谐 波 时 的 快速 
性 和 精度 优势 。 


2.0 


一 一 暂 降 幅 值 比 
i 电压 信号 (无 谐 波 ) 


0.02 0.07 0.12 人 0.22 0.27 0.32 
us 


(a) 暂 降 幅 值 比 


一 一 相信 天 赤 用 
a 电压 信号 (无 谐 波 ) 


0.02 0.07 0.12 OF 0.22 0.27 0.32 
ts 


(b) 相位 跳 变 角 
4 无 谐 波 干扰 条 件 下 的 二 点 法 电压 暂 降 检测 结果 
Fig.4 The results of voltage sag detection based on two 


points method without the harmonic interference 


加 入 谐 波 干扰 后 的 二 点 法 暂 降 特征 量 检测 结果 
如 图 5 所 示 。 可 见 ， 在 暂 降 信号 含有 谐 波 条 件 下 ， 
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暂 降 幅 值 比 及 相位 跳 变 角 检测 结果 存在 严重 误差 ， 
最 大 误差 分 别 可 达 42.25% 和 83.71%。 这 表明 二 点 
法 在 谐 波 干扰 环境 下 不 应 单独 使 用 。 


2.0 


一 一 暂 降 幅 值 比 
二 电压 信号 ( 含 谐 波 ) 


0 : 和 i nk |_350 
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ts 
(a) 暂 降 幅 值 比 
0.7854 T T 1 550 

一 一 相位 跳 变 角 9 
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(b) 相位 跳 变 角 
5 有 谐 波 干扰 条 件 下 的 二 点 法 电压 暂 降 检测 结果 
Fig.5 Theresults of voltage sag detection based on two 


points method with the harmonic interference 


4.4 滑 窗 小 波 与 二 点 法 电压 暂 降 检测 的 仿真 分 析 

最 后 ， 利 用 请 窗 小 波 与 二 点 法 对 仿真 信号 的 电 
压 暂 降 过 程 进行 检测 。 图 6 给 出 了 电压 暂 降 幅 值 比 
及 相位 跳 变 角 仿真 结果 。 为 便于 分 析 ， 仿 真 结果 图 
中 还 给 出 了 清 窗 小 波 得 到 的 电压 基 波 分 量 。 对 图 6 
分 析 可 知 ， 当 电压 发 生 暂 降 或 恢复 之 后 ， 一 个 工 频 
周期 之 内 ， 算 法 即 可 准确 地 跟 踩 暂 降 幅 值 比 及 相位 
跳 变 角 。 这 表明 本 文 方法 能 够 消除 谐 波 对 二 点 法 精 
度 的 影响 ， 提 高 暂 降 特征 量 计算 精度 。 尽 管 请 窗 小 
波 的 加 入 使 算法 产生 小 于 一 个 周期 的 延迟 ， 但 检测 
结果 已 可 满足 大 部 分 实际 工程 的 需要 。 
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一 一 牺 作 后 全 比 下 
区 滑 窗 小 波 得 到 的 电压 基 没 


0 1 1 1 1 ] i 
0.02 0.07 0.12 CE 0.22 0:27 0.32 
ts 


(a) 暂 降 幅 值 比 


201903.00116v1 


chinaXiv 


0.7854 T 


一 一 机 位 玩 变 角 9 
ee 滑 窗 小 滤 得 到 的 电压 基 波 


0.1963): 


0.3927F 


grad 


-0.9817| 


-1.5708 ， ， ， 
002 007 02 017 022 027 032 [3] 
LA 


(b) 相位 跳 变 角 
图 6 基于 滑 窗 小 波 与 二 点 法 的 电压 暂 降 检 测 结果 
Fig.6 The results of voltage sag detection based on sliding 


window wavelet and two points method 


表 2 给 出 了 不 同和 条 件 下 ， 基 于 背 窗 小 波 
与 二 点 法 计算 得 到 的 电压 暂 降 特征 量 平均 相对 误差 。 。 [4] 
分 析 可 见 ， 本 文 所 提 算 法 在 4 和 不 同 取 值 情况 下 ， 
均 能 得 到 高 精度 检测 结果 ， 不 存在 因 小 波 系数 最 大 
模 难 以 辨识 而 导致 的 暂 降 特 征 量 无 法 计算 问题 , 相 0] 
对 传统 小 波 变换 方法 而 言 ， 可 靠 性 优势 明显 。 

表 2 不 同和 上 + 条件 下 ， 基 于 滑 窗 小 波 与 二 点 法 的 

电压 暂 降 检测 误差 


Tab.2 The errors of voltage sag detection based on sliding 


window wavelet and two points method with the various [6] 


tiandt. 


序号 。 暂 降 时 刻 is ”恢复 时 刻 t/s ”误差 (%) 误差 (%) 


Casel 0.0766 0.2263 0.18 0.12 

Case 2 0.0781 0.2278 0.24 0.14 

Case 3 0.0797 0.2294 0.16 0.23 

Case 4 0.0813 0.2309 0.27 0.12 

Case5 0.0825 0.2322 0.24 0.25 中 
5 结论 


针对 小 波 变换 电压 暂 降 检测 算法 的 不 足 ， 提 出 
了 一 种 基于 清 窗 小 波 和 二 点 法 的 电压 暂 降 算法 。 实 
例 仿 真 表明 ， 所 提 算 法 不 存在 因 小 波 系数 最 大 模 难 
以 辨识 而 导致 的 暂 降 特征 量 无 法 计算 问题 ， 能 够 在 [8] 
谐 波 干扰 条 件 下 准确 地 计算 电压 暂 降 幅 值 比 及 相位 
跳 变 角 ， 检 测 结 果 延 迟 小 于 一 个 工 频 周期 ， 满 足 实 
际 工程 需要 。 
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